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In many biological systems, microorganisms swim through complex polymeric fluids, and usually 
deform the medium at a rate faster than the inverse fluid relaxation time. We address the basic 
properties of such life at high Deborah number analytically by considering the small-amplitude 
swimming of a body in an arbitrary complex fluid. Using asymptotic analysis and differential 
geometry, we show that for a given swimming gait, the time- averaged leading-order swimming 
kinematics of the body can be expressed as an integral equation on the solution to a series of simpler 
Newtonian problems. We then use our results to demonstrate that Purcell's scallop theorem, which 
states that time- reversible body motion cannot be used for locomotion in a Newtonian fluid, breaks 
down in polymeric fluid environments. 

PACS numbers: 47.63.Gd,47.63.mf,47.57.-s 

I. INTRODUCTION 

The physics of cell locomotion in viscous fluids affects many important biological processes [1], such as the journey 
of spermatozoa through the mammalian female reproductive tract [2], the mechanisms by which motile bacteria are 
able to progress towards high nutrient concentration [3], and the availability of plankton as food source for higher 
organisms in the ocean [4]. 

In many relevant instances, cells have to move through complex fluids, in particular during reproduction. In order 
to reach the uterus of the female and continue their journey towards the ovum, mammalian spermatozoa cells have to 
progress through the cervical mucus, a highly viscous and highly clastic cross- linked polymeric gel [2]. The rheology of 
cervical mucus depends on its hydration [5] , and varies during the female menstrual cycle [6] , but its typical viscosity 
is two to four orders magnitude larger than that of water [7-10], and its typical relaxation time, A, is in the 1 — 10 s 
range [5, 7, 10]. Since spermatozoa actuate their flagella with typical frequencies oj ~ 20 — 50 Hz [11], cell locomotion 
through the cervical mucus occurs therefore at high Deborah number, De = Xuj ^> 1, and elastic effects are expected 
to play a crucial role. 

Building on twenty years of research on the mechanics of locomotion in simple (Newtonian) fluids, Purcell detailed 
in his 1977 classical paper the physical principles of life at low Reynolds number [12]. In contrast, the basic properties 
of life at high Deborah number are not understood. Calculating the swimming speed of a given organism in a given 
complex fluid has only been solved for infinite models [13, 14], and the most basic questions remain unanswered: 
How different are the locomotion kinematics from those obtained in a Newtonian fluid? Can the nonlinear rheological 
properties of the fluid (in particular shear-thinning viscosity and normal stress differences [15]) be exploited to design 
new propulsion methods? 

Here, we address the problem of locomotion at high Deborah number analytically. We show that for small-amplitude 
swimming of a body in an arbitrary complex fluid, the swimming kinematics can be expressed as an integral equation 
on the solution to a series of simpler problems (motion in a Newtonian fluid) , thereby bypassing the explicit solution 
for the complete flow field. We then exploit our results to demonstrate explicitly that Purcell's scallop theorem — 
which states that time-reversible body motion cannot be used for locomotion in a Newtonian fluid [12] — breaks down 
in a polymeric fluid. 



* Electronic mail: elaugaOucsd.edu 



2 




FIG. 1: General statement of the swimming problem in a fluid: A body of fixed volume deforms its shape S(t) in a time-periodic 
fashion around an average shape, So- The surface deformation is prescribed in the swimming frame (Eulerian velocity u s ), and 
the unknown solid-body swimming kinematics (velocity, U; rotation rate, SI) are determined using the constraint of force-free 
and torque-free motion. 



II. NEWTONIAN SWIMMING 



We first recall the solution to the swimming problem in a Newtonian flow [16]. Consider an isolated three- 
dimensional swimmer of instantaneous surface S with normal n into the fluid. Lorentz' reciprocal theorem [17] 
states that for two arbitrary solutions of Newtonian Stokes flows with the same viscosity, (u, er) and (u, er), we have 
the equality 



u -&-iLdS= / / u-CT-ndS, (1) 

s J Js 

where u(u) and er(cr) are the velocity and stress fields. For (u, er) we consider the swimming problem (Fig. 1): In the 
swimming frame, the body prescribes its instantaneous surface velocity, u s , and as a result moves with instantaneous 
(but unknown) swimming velocity U and rotation rate SI, so that the surface velocity is given in the lab frame by 
u = U + f2 xx s + u ,s , for any point x s on its surface (in this paper, rotation rates and torques will be defined with 
respect to some arbitrary origin). For (u, er), we consider solid body motion of S with velocity TJ and rotation rate Ct, 
so that u = U + SI x x s on the surface. The body in the hat problem is therefore subject to an instantaneous force, 
F = jj a ■ ndS, and torque, L = jj x x (er • n) dS. Exploiting the fact that locomotion at low Reynolds numbers is 
force-free and torque-free, i.e. 



er • ndS = / / x x (er ■ n) dS = 0, (2) 

s J Js 

Eq. (1) leads to an equation for U and SI as 

F • U + L • SI = - J J n • e> • u s dS. (3) 

Eq. (3) states that, for a given shape (<S, n), and a given swimming gait (u s ), all six components of the swimming 
kinematics, (U, S~2) can be calculated using solely information about the dual problem of solid-body motion (F and 
L in Eq. (3) are arbitrary). Importantly, we note that the value of the fluid viscosity is irrelevant: As all hat terms 
in Eq. (3) are proportional to the viscosity, the relationship between the swimming gait (u s ) and the swimming 
kinematics (U, $7) is independent of the viscosity. 



III. LOCOMOTION IN NON-NEWTONIAN FLUIDS 



We now consider the case where swimming occurs in a complex fluid. The stress tensor, er, includes an isotropic part 
(the pressure, p), and a deviatoric component, r = cr+pl. We assume the velocity field, u, to be incompressible, and 
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therefore the equations for mechanical equilibrium in the absence of inertia are written as Vp = V • r and V • u = 0. 
For constitutive modeling, we assume that r can be written as a sum of different modes, r = J2i T ' S where each stress 
r l satisfies a nonlinear differential constitutive relationship of the form 

(l+A)T i + M i (T i ,u) = » ?i (l + iBi)7 + N i (7,u). (4) 

In Eq. (4), 7 = Vu + Vu T is the shear rate tensor, A4 and Bi are two sequences of linear differential operators 
in time representing polymer relaxation and retardation respectively, Mj and are two sequences of symmetric 
nonlinear operators representing transport and stretching of the polymeric microstructure by the flow, and r\i is the 
zero-shear rate viscosity of the i-th mode. The relationship between stresses and strain rates described by Eq. (4) is 
a very general differential constitutive relationship [15, 18-21], which includes as particular cases all classical models 
of polymeric fluids [28]. 

We consider a body performing periodic small-amplitude swimming motion in a fluid described by Eq. (4). Its 
undcformed surface shape is termed So, parameterized by xjj , and we define e as the amplitude of the periodic surface 
distortion non-dimensionalizcd by a typical swimmer length (e <C 1). Material points on the swimmer shape, x s , 
are assumed to display time- variations of the form x. s (x^,t) = Xg + exf (xq,£), and the function xf is assumed to 
be periodic in time with period T. Such Lagrangian boundary motion forces the fluid to move through the no-slip 
boundary condition, u s (x s ) = dx. s jdt. 

We solve the swimming problem as a domain perturbation expansion, where the fields of interest are written as 
regular perturbation expansions, with boundary conditions expressed on So- Specifically, we write 

{u,T,p,er} = e{ui,ri,pi,cr 1 } + e 2 {u2,-r 2 ,_P2,CT2} + ••• (5) 

which are all functions of (x, t) , and are defined on the zeroth-order surface So ■ The boundary condition for the 
surface velocity reads 

u s = euf(x 9 ,t) + e 2 uf(x 9 ,t) + ... (6) 
and swimming occurs with the kinematics 

{U,0} = e{U 1 ,ft 1 } + e 2 {U 2 ,n 2 } + ... (7) 

so that on the swimmer surface we have u„ = U„ + tl n x xjf + uf , for n = 1,2, .... Based on the Newtonian case, we 
expect to obtain no swimming at order e, but non-zero time-averaged locomotion at order e 2 [22]. 

A. First-order solution 

At order e, the constitutive model, Eq. (4), is linearized 

(1+A)ri =%(1+Si)7i, (8) 

associated with boundary conditions uf = dxf/dt, evaluated at (xg,i). Since the surface motion is time-periodic, 
we introduce Fourier series, and write, for any field f(t), f(t) = X^oo f^e mojt where u = 2n/T and = 
T So f{t)e- mut dt. In Fourier space, Eq. (8) then becomes 

fi' ( " ) (x) = ft(n) 7 ^(x), (9) 

where Qi(n) is the i th relaxation modulus of the n th Fourier mode. Since we have r = T ' \ we § e t the constitutive 
equation for the total first order deviatoric stress as 

^(x) = S(n) 7i n) (x), Q(n) = $>(n). (10) 
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We see from Eq. (10) that, for each Fourier mode, the swimming problem is a Newtonian problem with a complex 
viscosity ((?). We have to solve Vp^ = Q(n)V 2 VL( l \ V • = 0, subject to the boundary condition u^(xjf) = 
+ fij™'' x Xq + uf ''^ (xq ). Applying Eq. (3), we obtain the swimming kinematics for each Fourier mode 

F-Ui n) +L-nS n) = - ff n .^.uf (n) (x s )d5. (11) 

J J So 

Since the value of the viscosity for the hat fields in Eq. (11) is arbitrary, we can take it to be some fixed reference 
viscosity. In addition, as So does not depend on time, we can Fourier-invert Eq. (11) to obtain the locomotion in the 
time domain 

F-Ui(i)+L-ni(t) = - // n -<T-uf(x^,t)dS. (12) 

J J So 

The solution at order e leads thus to the same swimming kinematics as in a Newtonian flow (Eq. 3). In addition, 
since uf = <9xf /dt, we get that (uf ) = 0, where (.) denotes time- averaging over one period of body deformation (i.e. 
the zeroth Fourier mode). From Eq. (12) we therefore sec that (Ui) = (fii) = 0. As in the Newtonian case, there 
is no time-averaged locomotion at leading order, and swimming is quadratic in the amplitude of the surface motion 
[22]. 



B. Second-order solution 

At order e 2 , the constitutive relationship for each mode, Eq. (4), is written as 

(1 + A t )r 2 = 7?,(1 + Bi)j 2 + Hi[ui], (13) 

with 

Hj = m ■ { 7l : [(V^V U N,)] - r\ : [(V^V U M,)]} (14) 

where the gradients in Eq. (14) are evaluated at (0,0), and with t\ and ■y 1 related through Eq. (8). Since we are 
interested in the time- averaged swimming motion, which we expect occurs at 0(e 2 ), we now consider only time- 
averaged quantities. Averaging Eq. (13) leads to 

(rj>=w(72> + (H i [u 1 ]) ) (15) 

and therefore the time-averaged stress is given by 

(<T2> = -<P2>l+r?(72> + <S[ui]), (16) 

where V^Y^iVi and ( s [ u i]> = Ei( H 4 u i]>- 

To derive the swimming kinematics, we apply the principle of virtual work using the following two problems: (i) 
solid body motion of the shape So in a Newtonian fluid of viscosity rj (the same viscosity as in Eq. 16) , with velocity 
and stress fields given by u and <r and (ii) time-averaged swimming with flow velocity (u 2 ) and stress field (er 2 ) given 
by Eq. (16). 

Since mechanical equilibrium is written V ■ <x = V ■ (er 2 ) = 0, we have equality of their dot products with the 
opposite velocity field, [V • a] ■ (u 2 ) = [V • (er 2 )] • u, and integration over the volume of fluid Vq outside of So leads to 

ff n - a- (u 2 ) dS - ff n • (<x 2 ) • udS (17) 

J J Sq J J So 

= Iff (<t 2 ): VudV- fff &:W(u 2 )dV, 
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where we have used integration by parts, and the fact that no is directed into the fluid. If we then insert Eq. (16) 
into the right-hand side of Eq. (17) wc obtain 

(18) 



<r 2 ) : VudV - fff a: V(u 2 ) dV 

J J Jv a 

fff (S[ui]):VudV, 
J J Jv 



and the Newtonian components of both <x and (<r 2 ) have disappeared due to symmetry and incompressibility. Con- 
sequently, Eq. (17) becomes 



n • <x • (u 2 ) dS ■ 

S[ui]) : \7udV, 



n • (er 2 ) • udS 



(19) 



and only the deviation from Newtonian behavior, X, remains in the integral formula. This result is reminiscent of 
past work quantifying small viscoelastic effects on particle motions [23]. 

On the surface So we have (u 2 ) = (U 2 ) + (0 2 ) x Xq + (uf ), where a Taylor expansion of the boundary conditions 
around xf leads to uf (xf,i) = — xf • Vui, so that Eq. (19) becomes 



F-<U 2 )+L-<n 2 ) 



n -6--(uf)dS* 



(20) 



So 



+ 



// n - (<r 2 ) -ud5+ fff <E[ui]> : VudV. 

J J So J J JVo 



The final step in the calculation consists in enforcing the force-free and torque-free condition for the swimmer. On 
So we have u = U + ft x xjf , so that 



/ / n • (er 2 ) • udS = 

J J So 



ff n -(<r 2 )d5 

J J So 



■U 



(21) 



+ 



So 



xjf X (n • (er 2 ))dS 



n. 



The terms in brackets in Eq. (21) are related to the second forces and torque on the swimmer at order e 2 , and 
can be evaluated using differential geometry. Let us write the time-varying shape of the swimmer as x s = xjf + 
eno#i(xo , t) + where the function Si, with units of length, represents the normal extent of the surface deformations. 
When Si = 0, the shape of the swimmer does not change with time, and all surface motion is tangential (u s • n = 0, 
so-called squirming motion), whereas for Si ^ the body also undergoes normal deformation and varies its shape 
periodically. If we write the normal to the surface as n = no + eni + differential geometry considerations leads to 
the evaluation of the force, F 2 , and torque, fl 2 , on the swimmer at order e 2 , as given by 



F 2 - 



So 



ni • <Ti + n • er 2 + Si 



fh = [[ x s 

J J So 



dcri 

dn 



ni • crx + n • er 2 + Si 



dS, 

dcri 

dn 



dS 



(22a) 
(22b) 



where d/dn = n • V denotes the normal derivative to So- Since locomotion occurs with no force or torque, we have 
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F 2 = f2 2 = 0, and therefore after taking time-averages of Eq. (22), we obtain 



// n -(<x 2 )dS 

J J So 



So 



(ni • <ti) + n • ( 5i 



l l So ° 

l l So ° 



x (n • (<r 2 )) dS 1 



<9<Xl 

dn 



dS, 



(ni • <7i) +n • ( 5i 



do - ! 

9n 



(23a) 



(23b) 



dS. 



C. Life at high Deborah number 

To obtain the final integral formula, we insert the result of Eq. (23) into Eqs. (20) and (21) to obtain the integral 
relationship 



F • (U 2 ) + L • (n 2 ) = 

-// no-(T-<uf)dS+ [[[ <£[m]> : VudV 

J J So JJJVo 



(24) 



(ni ■ <Ti) + n • ( S 



Vo 

dcri 

dn 



dsy-u 



Xn X 



So 



(ni • o-i) + n ■ ( Si 



dcri 

dn 



The result expressed by Eq. (24) is the non-Newtonian equivalent of the Newtonian integral formula, Eq. (3). It shows 
that one can compute the time-averaged swimming kinematics for locomotion in a complex fluid, using knowledge of 
a series of simpler problems. Indeed, to compute (U 2 ) and (f2 2 ) from Eq. (24), and beyond the necessary knowledge 
of the surface motion of the swimmer, one needs to know the velocity and stress held for solid body motion of So (i.e. 
the fields u and <x), and the velocity and stress held for the hrst-order solution (i.e. Ui and cri). As discussed above, 
and shown in Eq. (10), the hrst order solution can be found in frequency space by solving a series of Newtonian flow 
problems. Consequently, the computational complexity to evaluate the terms in Eq. (24) is that of a succession of 
Newtonian flow problems, and therefore using this method one bypasses entirely the calculation of the second-order 
flow and stress field. Notably, the hnal result can be applied to flows with arbitrary large Deborah numbers, as is 
relevant in cell locomotion. Note also that for squirming motion of the sphere, for which Si = and ni = 0, Eq. (24) 
is greatly simplified. 



IV. BREAKDOWN OF THE SCALLOP THEOREM 

As an application of our results, we demonstrate that Purcell's scallop theorem [12] breaks down in a polymeric 
fluid. We consider the axisymmetric squirming motion of a sphere (radius, a) in an Oldroyd-B fluid [15, 18-21]. 
Purcell's scallop theorem states that if the surface motion is time-reversible, we have (U) = (ft) = and therefore 
the Newtonian contribution to Eq. (24) averages to zero, 

/ / n • & ■ (uf ) dS = 0. (25) 

J J So 

In addition, we consider axisymmetric surface deformation so that we have (ft 2 ) =0. As a consequence, the integral 
equation leading the average swimming speed, Eq. (24), simplihcs to 

F-(U 2 )= /// <E[ui]):VudV. (26) 

J J JVo 
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FIG. 2: Distribution of surface velocity, vm, (arbitrary units) for axisymmetric spherical squirming motion and breakdown of 
the scallop theorem, (a): Surface velocity as a function of the polar angle, 6; (b): Color map of the the surface velocity. Both 
figures illustrate the fore-aft asymmetry of the tangential surface motion. 



For constitutive modeling, we consider an Oldroyd-B fluid, which represents a polymeric fluid as a dilute solution of 
elastic dumbbells [21], and for which the relationship between stresses and rate of strains is given by 



t + Ai r= rj[j + A 2 7], 



(27) 



where a= da/dt + u • Va — (Vu T • a + a • Vu) is the upper-convected derivative for the tensor a. In Eq. (27), Ai and 
A2 are, respectively, the relaxation and retardation time scales for the fluid. If rj s denotes the solvent viscosity, and 77 
the polymer viscosity, we have A2/A1 = r\ s jr\ < 1. 

For a time-reversible deformation, we consider a simple sinusoidal gait of the form uf (xg,t) = vj?(xQ ) cos cut, so 



that u 



S,(n) 



(xo S ) 



for all n ^ ±1, and uf (±1) (xj?) 



v]?(xq)/2 otherwise. At order e, only the Fourier modes with 



n ^ ±1 are non-zero, and we have 



= S*7i 



G = v 



1 + i\2UJ 



(28) 



1 + iXiut ' 

where {.}* denotes the complex conjugate. The spatial distribution of surface deformation, described by vn , is assumed 
to be axisymmetric. The symmetry axis is denoted e 2 (see Fig. 2), the polar angle is 9, and is associated with the 
orientation vector eg. In the frame moving with the swimmer, we prescribe 

vf(xf) =3awsm0(l + cos0)ee, (29) 

which is illustrated in Fig. 2. Note that the velocity distribution described by Eq. (29) is fore-aft asymmetric, which 
is necessary in order to obtain net locomotion with an actuation varying sinusoidally in time. 

Given Eq. (29), we can then calculate the unsteady swimming at order e from Eq. (12), and we find Ui = 
2au> cosute z . As a result, the surface distribution of velocity in the lab frame is given by Ui(xq,£) = V| (xq ) cos uit, 
where 



II I 



2auj cos 9 e r + auj sin 9(1 + 3 cos 9) eg 



(30) 



Given Eq. (11), it is then easy to show that each Fourier component of the entire flow field is identical to that obtained 
in the Newtonian problem. Consequently, if V||(x) denotes the Newtonian velocity field associated with the lab-frame 
boundary conditions V||(xq) on So, we obtain at first order Ui(x, t) = V|| (x) coswt. The velocity field V||(x) with 
boundary conditions from Eq. (30) can be found using the Legendre polynomials method pioneered by Blake [24], 
and we get V|| = V||, r e r + v\\ t geg with 



'II,- 



v \\,0 



a 

2— cos( 

„3 



(3 cos^-l) --- 



o^ 4 ■ 
3— sm( 

r 4 



cos( 



(31a) 
(31b) 
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At order e 2 , straightforward algebra allows us to obtain the deviation from Newtonian behavior, in Eq. (16), as [13] 



<s[ Ul ]) = x 



2(1 + De 2 ) 
v r V7||-(Vv[-7||+7 r Vv||) 



(32) 



where De = Xilj is the Deborah number for the flow. 

Finally, the hat problem in Eq. (26) is the solid body translation of the sphere, with velocity field given by [25] 



. 3 
u = -a 

4 



• U + la 3 

4 



3rr 



U 



(33) 



together with Stokes law, F = — 67r?7oU. 

By symmetry, we expect that average swimming will occur along the z direction, so that (U2) = (U2)e z and by 
choosing U = Ue Z} the left-hand side of Eq. (26) is given by —Qir-qaV (U2) ■ Given Eqs. (32), (31) and (33), we can 
evaluate the right hand side of Eq. (26) and obtain 



fff 1Tr 2 2rrV(^i - A 2 ) 299tt 

/// (S[ui]) : Vudy = a 2 u; 2 U 1 

J J JVn 



(34) 



Recalling that A 2 = 
squirming sphere as 



ivo 1 + Bc z 25 

Xir/s/r], we obtain the explicit formula for the time-averaged swimming speed, (U2), of the 



(U2) = 



De 



1 + De 2 



-1 A 



V 



(35) 



where A = 299/150 w 1.993. The result of Eq. (35) demonstrates explicitly that the scallop theorem breaks down in 
an Oldroyd-B fluid: The swimming gait is a sinusoidal function, and therefore time-reversible, yet the force-free body 
swims on average. In the Newtonian limit where De = 0, we have (U2) = and the result of the scallop theorem is 
recovered. Note that since rj s < rj, we have {U2) < 0. High surface shear is localized on the top of the sphere (see 
Fig. 2b), so this is also where high normal-stresses differences are localized, and the sphere is being pushed from the 
top to swim in the — z direction. 



V. PERSPECTIVE 

In this paper, we have addressed the most basic problem in the locomotion of microorganisms: For a given swimming 
gait, at which speed is the organism expected to swim? The solution to this problem is known in the case where 
the fluid is Newtonian, and given by Eq. (3), but is not known for complex polymeric fluids displaying a nonlinear 
relationship between stress and strain rates. We have considered the time-periodic small-amplitude locomotion of a 
dcformable body in an arbitrary complex fluid. We have shown that the time-averaged swimming kinematics of the 
body (translation and rotation) are given by an integral formula on a series of simpler Newtonian problems. The 
final formula, Eq. (24), can be applied for high Deborah numbers, which is the relevant limit for the locomotion of 
swimming cells in mucus, and provides the first formal framework to address locomotion in complex fluids. In addition, 
our results are valid beyond the biological realm, and can be used in particular to the quantify the locomotion of 
synthetic micro-swimmers [26]. 

As an application of our results we have constructed an explicit example of a deformable body that swims using 
a time-reversible stroke in a polymeric fluid. This example demonstrates formally the breakdown of Purcell's scallop 
theorem in complex fluids for a finite-size, force-free and torque-free swimmer. Note that the final formula for the 
time-averaged swimming speed of the body, Eq. (35), is reminiscent of recent work on the force generated by flapping 
motion in polymeric fluids [27]. The implication of this result, more generally, is that it is possible to exploit nonlinear 
rheological mechanisms (in our case, the existence of normal-stress differences) to design new swimming methods. 
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Finally, we note that recent work on infinite models for swimmers deforming in a wave-like fashion showed that, 
for a given swimming gait, swimming is always slower in a polymeric fluid than in the Newtonian limit [13, 14]. The 
final integral formula for the swimming speed we obtain here, Eq. (24), explicitly shows that in general the beneficial 
vs. detrimental impact of the polymeric stresses on the swimming performance cannot be established a priori. 

The results above could be extended in many different ways. In particular, the method of expansion outlined in 
the paper could be further continued, and all Fourier components of the flow at higher order in the amplitude of 
the surface deformation could be formally calculated. Similar work could also be performed near boundaries, or in 
the presence of other swimmers, and therefore could be exploited to characterize the effect of polymeric stresses on 
collective locomotion. The application of our results to different swimmer geometries and various modes of surface 
swimming, including flagella-based, will be reported in future work. 
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